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A Fully Many-Body Localized (FMBL) quantum disordered system is characterized by the emer¬ 
gence of an extensive number of local conserved operators that prevents the relaxation towards 
thermal equilibrium. These local conserved operators can be seen as the building blocks of the 
whole set of eigenstates. In this paper, we propose to construct them explicitly via some block 
real-space renormalization. The principle is that each RG step diagonalizes the smallest remaining 
blocks and produces a conserved operator for each block. The final output for a chain of N spins is 
a hierarchical organization of the N conserved operators with layers. The system-size nature 

of the conserved operators of the top layers is necessary to describe the possible long-ranged order 
of the excited eigenstates and the possible critical points between different FMBL phases. We dis¬ 
cuss the similarities and the differences with the Strong Disorder RSRG-X method that generates 
the whole set of the 2^ eigenstates via a binary tree of N layers. The approach is applied to the 
Long-Ranged Quantum Spin-Glass Ising model, where the constructed excited eigenstates are found 
to be exactly like ground states in another disorder realization, so that they can be either in the 
paramagnetic phase, in the spin-glass phase or critical. 


I. INTRODUCTION 

The field of Many-Body Localization (MBL) has attracted a lot of interest recently (see the recent reviews 0,0 and 
references therein). The general goal is to understand the unitary dynamics of isolated random interacting quantum 
systems and to determine how they can avoid thermalization and remain non-ergodic. 

One of the most important characterization of the Many-Body localized phase is that excited eigenstates display an 
area-law entanglement 0 instead of the volume-law entanglement of thermalized eigenstates. This property has been 
used numerically to identify the MBL phase Si and to show the consistency with other criteria of MBL 0. The 
fact that excited eigenstates in the MBL phase are similar to the ground-states from the point of view of entanglement 
suggests that various approaches that have been developed for ground states in the past can be actually adapted to 
study the MBL phase. One first example is the efficient representation via Density-Matrix-RG or Matrix Product 
States and Tensor Networks [l0. Another example is the Strong Disorder Real-Space RG approach (see pm 
for reviews) developed by Ma-Dasgupta-Hu [l0 and Daniel Fisher IT5l| to construct the ground states of random 
quantum spin chains, with its extension in higher dimensions d = 2, 3,4 |16l426j| . This approach has been extended 
into the Strong Disorder RG procedure for the unitary dynamics [^ . , and into the RSRG-X procedure in order 

to construct the whole set of excited eigenstates [27l - P . It should be stressed that these two Strong Disorder RG 
procedures based on the spin variables are limited to the MBL phase, whereas the current RG descriptions of the 
MBL transition towards delocalization are based on RG rules for the entanglement or for the resonances [s^. 

Another recent essential idea in the field is the claim that a fully many-body localized system (FMBL) can be 
characterized by an extensive number of emergent localized conserved operators [36l - l44| . More precisely for a random 
quantum chain containing N spins ai described by Pauli matrices, it should be possible to introduce N localized 
pseudo-spin operators rf that commute with each other and with the Hamiltonian. In terms of these localized 
pseudo-spins r^, the Hamiltonian can be rewritten as 

N k 

jjfmbl ^ ( 1 ) 

k—{) l<i\<i 2 <..<ik<N q—1 

containing 2^ independent couplings . It should be stressed that the diagonalization of any Hamiltonian with 

2^ energy levels can be rewritten as Eq. [T] since the 2^ couplings are sufficient to reproduce the 2^ energies. So the 
non-trivial statement in the FMBL Hamiltonian of Eq. [T] is that the pseudo-spins rf are obtained from the real spins 
uf by a quasi-local unitary transformation, and that the couplings J* decay exponentially with the distance with 
a sufficient rate (see [36l - l4l| and the reviews 0i for more details). On the contrary, in the delocalized phase, the 
pseudo-spins t? that diagonalize the Hamiltonian are delocalized over the whole chain. 

Since the essential property of these conserved pseudo-spins rf in the FMBL phase is their locality, it seems natural 
to try to generate them locally and homogeneously in space. In the present paper, we thus propose to construct them 
via some block Real-Space Renormalization procedures that generalize previous procedures introduced for ground 
states. The principle is that at each RG step, instead of projecting always onto the lowest-energy subspace, one 
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keeps the projection onto the lowest energy subspace and the projection onto the highest energy subspace, so that the 
choice between the two subspaces corresponds to a conserved pseudo spin. In spirit, the idea is thus very close to the 
RSRG-X method where each decimation produces a bifurcation in the spectrum and identifies a conserved 

pseudo-spin. In practice, the output is however somewhat different. For instance for a chain of N spins with 2^ 
eigenstates : 

(i) the RSRG-X method produces a binary branching tree with N layers, corresponding to the N successive 
decimations, so that at the end of the procedure, the 2 ^ leaves of the tree are the eigenstates. 

(ii) the present block Real-Space approach produces a hierarchical tree with layers for the conserved operators 

themselves : the first RG step produces in parallel ^ conserved operators associated to the blocks of L = 2 sites, the 
second RG step produces ^ conserved operators associated to the blocks of L = 4 sites, and so on, i.e. the k RG step 
produces ^ conserved operators associated to the blocks of L = 2^ sites. So the first RG steps produce an extensive 
number of very local conserved operators as expected. However the few conserved operators associated to the last 
RG steps are system-size and not ’local’ anymore. But their system-size nature is necessary to describe the possible 
long-ranged order of the excited eigenstates and the possible critical points between different FMBL phases . 

The paper is organized as follows. We consider the random quantum Ising chain with random transverse fields and 
with random couplings that are either short-ranged in section m or long-ranged in section Hill Our conclusions are 
summarized in section IlYl 


II. RANDOM QUANTUM ISING CHAIN 


The quantum Ising chain 


H = - E (2) 

i i 

with random transverse fields hi and random nearest neighbor couplings Jiy+i is ’trivially’ in the Fully Many-Body 
Localized phase, as a consequence of the Jordan-Wigner mapping onto free fermions. Nevertheless, we feel that it 
is useful as a toy model to see how the present approach works in the simplest possible case. In this section, we 
thus describe how the Fernandez-Pacheco self-dual block-RG used to construct first the ground state of the pure 
chain [i^, and then the ground state of the random chain [H-Elj can be extended to construct the whole set of 
eigenstates. 


A. Diagonalization of the intra-block Hamiltonian 


For each block of two spins {2i — 1, 2z), the intra-block Hamiltonian chosen by Fernandez-Pacheco [4 


H, 


(2i-l,2i) 


— —h2i-lO'2i-l — J2i-l,2i0'2 




(3) 


commutes with cr^j. So for each eigenvalue 821 = ±1, the diagonalization of the effective Hamiltonian for the remaining 
odd spin a 2 i-i 


h2ili{S2i) — —h2i-l(J2i-i — J2i-l,2iS2iCr2i-l ( 4 ) 

leads to the two eigenvalues (that are independent of the value 821 = ±1 as it should by symmetry) 

'^2i-l('^2i) = 2i (3) 

with the following corresponding eigenvectors 


|A2-_i(^2*)> = 




1 + 


J 2 i-lS 2 i 




■|'S'2i-l — +l,<S'2i > -F^ 


1 - 


J 2 i-lS 2 i 




■ |5'2i-l — —1, 82 i > 


■\ 


1 - 


^ 2 i-lS 2 i 


-\82i-1 — 


Id- 


•J2i-lS2i 




-\82i-1 — — 




2 


+ 1 ) 821 ^ 4 "^ 


2 


1 ,^ 2 . > ( 6 ) 
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To construct the ground state one chooses to keep only the two degenerate lowest states |A {S 2 i = ±) >. 

Here, as in the RSRG-X method 2j31|, we wish to keep also the two degenerate highest levels |A+(5'2i = ±) >. 
It is thus convenient to introduce a pseudo-spin for the choice between these two possible energy levels and a 

spin (J^( 2 i) to distinguish the two degenerate states within each energy level. More precisely, the four eigenstates are 
relabeled as follows 


^2i-l( + ) > 

= kfl(2i) 

= +: ^R{2i) = + > 


'^2z-l(~) > 

— kfl(2i) 

= U'fl(2i) = “ > 


|A2Vi(-) > 

= l'^fl(2i) 

= ®’fl(2z) = ~ > 


A2i-l(+) > 

— kfl(2i) 

= ^R{2i) = + > 

(7) 


so that the intra-block Hamiltonian of Eq. [3] depends only on the pseudo spin T/j( 2 i) 


with the field 


H(2i-l,2i) — -^Ri2i)T^{2i) 


( 8 ) 


n 


R{2i) 




'^2i-l,2i 


(9) 


Since this procedure is applied independently to all blocks (2i — l,2i), the total intra-block Hamiltonian reads in 
this new basis 


j^intra 


= E^ 


intra 

{2i-l,2i) 


i 


= - Vh 


7 ^‘‘R(2i)T}i{2%') 


( 10 ) 


For a chain of N spins, this intra-block Hamiltonian depends only on the pseudo-spins whose flippings are 

associated to the ’large’ fields ^R( 2 i), but is completely independent of the ^ spins (jR( 2 i), whose flippings have ’no 
cost’ yet when one takes into account only In the following, we thus make the approximation that the pseudo¬ 
spins T^( 2 i) * = li2,Y are the first ^ emergent local conserved operators. This approximation is possible only 

when the initial transverse fields hi and/or the initial couplings J^y+i are distributed with continuous probability 
distributions, so that there isn’t any exact degeneracy between the fields H/j( 2 i) of Eq. (For instance for the pure 
chain where all fields ^R[ 2 i) coincide, it is clear that this approximation would be meaningless). 


B. Renormalization of the extra-block Hamiltonian 


The extra-block Hamiltonian defined in terms of the initial spins 


jje.tra = i/ _ ^ J2,_2<_2a2Vl (H) 

i i 

has to be rewritten in terms of the new spin variables (T^( 2 i) > o'R( 2 i) ) introduced in Eq. [T] Within the approximation 
that the pseudo-spins T^( 2 i) emergent local conserved operators, one neglects all terms involving the flip operators 
obtains the following effective Hamiltonian for the remaining spins cr^( 2 i) 


h: 


eff 




R(‘ 2 i)^R{ 2 i) 


E '^fl(2i-2),tj(2i) )^R(2i-2) ^R(2i) 

i 


( 12 ) 


where the renormalized transverse field 


hR(2z) - h2i 


\j^2i-l + '^2i-l,2 


(13) 


has decreased h^( 2 i) < h 2 i and where the renormalized coupling 


jR{2i-2),R(2i)iT^{2i)) “ TR(^2i)J'2i-2,2i-l 


J2i-l,2i 




2i 


(14) 


has also decreased in absolute value \JR( 2 i- 2 ),R{ 2 i)\ < \J 2 i- 2 , 2 i-i\- The renormalized fields hR(^ 2 i) uud couplings 
JR(2i-2),R(2i) are thus typically smaller than the initial fields and couplings, and thus also typically smaller than the 
fields ^R{ 2 i) of Eq. |9] associated to the conserved pseudo-spins T^( 2 i) ■ 

Eqs [13] and [T3| are the same as in the RG rules for the ground state except for the presence of the pseudo 

spin T^( 2 i) change the sign of the effective coupling between the spins o-R(^ 2 i- 2 ) ^R( 2 i)- 
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C. Comparison with the RSRG-X method 

The differences with the RSRG-X method are the following : 

(i) here one makes an a-priori arbitrary choice with the couplings J 2 i-i, 2 i) concerning the odd spins, while 

the RSRG-X method selects the biggest coupling of the whole chain at each step. 

(ii) here the renormalized field n_R( 2 i) = sj 2 i appear also in the eigenstates of Eq. [51 and in 

the RG rules of Eos [T51 and UTl take into account any ratio h 2 i-i/J 2 i-i, 2 i between the transverse field /i 2 i-i and the 
coupling J 2 i-l, 2 i^ whereas in the RSRG-X method it is approximated by the maximal coupling, namely either /i 2 i-i 

or J2i-l,2i- 

As for the RSRG method for the gound state 0 [ll, the RSRG-X method [13 is expected to become asymptotically 
exact at large RG steps when the RG flow is towards an Infinite Disorder Fixed Point, which happens at the critical 
point between two many-body-localized phase [l^- However for the present goal to construct the extensive local 
conserved operators of the MBL phase, it is clear that the majority of the pseudo-spins are produced by the first 
RG steps on small scales, and not by the asymptotic RG flow at large scales. In a previous work concerning the 
Shannon and Renyi entropies of the ground state, that are also dominated by the first RG steps on small scales, it 
was found that the block RG approach gives a better approximation of the multifractal dimensions than the RSRG. 
So here also we expect that the better RG rules of (ii) overcompensate the drawback of the arbitrary choice of (i) and 
that the block RG will produce a better approximation for the whole set of local conserved operators. 


D. Reration of the renormalization process 

In summary, at the end of the first RG step, the spectrum of the 2^ levels of the chain containing N spins has 
been decomposed into 2^ groups of 2^ levels each. A given group is labeled by the values of the ^ pseudo-spins 
T^( 2 j) = ±1 and has for reference energy the contribution of Eq. fTOl 

271-1 

i=l 


N 

With respect to this reference energy, the 2 2 levels are those of the quantum Ising chain described by the effective 
Hamiltonian of Eq. [T^] for the ^ spins with the smaller renormalized fields /i_R( 2 i) < and the smaller 

renormalized couplings | Jfl( 2 i- 2 ),fl( 2 i)I < |-^ 2 i- 2 , 2 i-i|- 

N N N 

The RG procedure may be iterated as follows. Each group of 2~ levels will be subdivided into 2~ sub-groups of 2~ 
levels each. Each subgroup will be labeled by the ^ pseudo-spins T^( 2 i) = il of the first generation that define the 

group and by the subsequent choice of the ^ pseudo-spins ''‘^ 2 ( 4 ^ = ±1 of the second RG step, so that the reference 
energy of this subgroup is now 


27.-1 2"-2 

ER'^i{'^R{2t)}'A'^R2{4i)}) = “ 51 “ 55 (4i) (4i) 

i=l i=l 


(16) 


where the fields of the second generation 


n 


R^(4i) 




■ 2 ) 


R{2i-2),R(2i) 


(17) 


are actually independent of the multiplicative sign appearing in the renormalized coupling in Eq. [TT] However 

the second generation pseudo-spin T^ 2 ( 4 j) itself depends on the two pseudo-spins values 'r^( 4 j_ 2 ) = and T^( 4 q = ±1 
of the first generation. Of course, if one insists on reproducing the form of Eq. [U it is always possible to rewrite any 
function of these two pseudo-spins as the polynomial 


fAR(4i- 


R(4i-2)^ ‘R(4i) 


)= J ( 0 . 0)+^(1 




' R{4i-2) 


- 2 ) + / 


( 0 . 1 ), 


^R(4i) 


+ f 


( 1 , 1 ) 


'’'R(4i-2)'’'R{4i) 


(18) 


where the four coefficients are chosen to reproduce the only four possibles values /(±1,±1). But in the present 
framework, it is simpler to stick to the pseudo spin T-^ 2 ( 4 q as defined by the renormalization process in each sector, 
and to keep in mind that it depends on the lower-generation pseudo-spins within its block. 
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__ ;v 

With respect to the reference energy of Eq. [161 the 2~ levels of this subgroup are those of the quantum Ising chain 
of ^ spins cr|j 2 (- 4 j) with the appropriate renormalized transverse field and coupling obtained by the next iteration of 
the RG rules of Eq. [13] and [14] 




J, 


T-R{4i-2) 


^K(4i-2) + ’^R{4i-2),R(4i) 


R^(4i-4)M^(4i) 


~ TR2(4i)’^Ri4i-4),R{4i-2) 


jR(4i-2),R{4i 


^R{4i-2) '^R(4i-2),R(4i) 


(19) 


For a chain of IV = 2" spins, the procedure ends after the RG step : the choice of 2"“^ pseudo spins T^( 2 i) 
the first generation has been followed by the choice of 2"“^ pseudo spins 'r^ 2 ( 4 j) of the second generation, then by the 
choice of 2"“^ pseudo spins T^ 2 (gj) of the third generation, etc... up to the choice of 2° = 1 pseudo spin of fhe 

last generation; there remains a single spin cr/{n( 2 '») with its renormalized field h/jTi( 2 »») so that the Hamiltonian is 
fully diagonalized and reads 


2’*“! 2”“^ 2"“^ 

HRn = - ^ ^R{2i)TR(2t) - X! “ X! ^iSi)'''R^(Si) 

i—1 i—1 i—1 

■■■ — f^fl"(2")'r^"(2n) ~ ^fl"(2'*)0'fl»(2") (20) 

i.e. o'|jn( 2 n) is the last conserved operator. 

Note that Eq. [501 contains only 2'^ = N parameters (namely (2” — 1) parameters of the various generations 
1 < k < n and one parameter h/jn( 2 -)) instead of the 2 ^ coefficients of Eq. [Tjto parametrize the 2 ^ energies : this is 
a consequence of the free-fermion nature of the quantum Ising chain of Eq. [21 For models that cannot be reduced to 
free-fermions, the diagonalized Hamiltonian will be more complicated than Eq. 1201 as shown explicitely in section iHll 


E. Properties of the eigenstates 


In terms of the ratio 




Ji—\A 


( 21 ) 


the RG rules of Eq. [T3| and [T4| reduce to the simple multiplicative rule 


K 


R(2i-2) = 


jR 

R{2i-2),R(2i) 

hR{ 2 i) 


- TR{2i)^‘2i-2K2i-l 


( 22 ) 


The conclusions for the amplitudes |iV| are thus exactly the same as for the ground state case corresponding to the 
choice T^( 2 i) = +1 discussed in detail in previous works : 

0 and the eigenstates are paramagnetic. 

^ -foo and the eigenstates display a long-ranged 


(i) in the region In hi > In the flow is attracted towards K 


(ii) in the region In hi < In Ji^i+i, the flow is attracted towards \K\ 
order adapted to the signs = ±1. In particular in the middle of the spectrum, as a consequence of the random 
signs = ±1, the eigenstates are in the spin-glass phase for any initial distribution of the couplings (even if the 
initial model is ferromagnetic Ji 4+1 > 0). 

(iii) at the phase transition In hi = In the eigenstates are critical and described by the Infinite Disorder Fixed 

point with the activated exponent ip = 1/2, with the typical correlation length exponent vum = 1 and the finite-size 
correlation exponent i/ps = 2, in agreement with the Fisher Strong Disorder exact results |15l |. 

(iv) the eigenstates have also exactly the same Shannon-Renyi entropies as the ground state studied in [^ . since 
the signs r/ completely disappear from the RG rules for the Shannon-Renyi entropies. 

The entanglement properties of the excited eigenstates are also exactly the same as for the ground state, namely 
an area law (i.e. a constant here in d = 1) outside criticality, and a lo gari thmic violation at criticality, whose exact 
behavior has been computed via the Strong Disorder RG approach ^.153 1. 
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III. LONG-RANGED DYSON QUANTUM ISING MODEL 

Many-Body Localization has been studied for various long-ranged power-law interactions . Any long-ranged 

model has a Dyson hierarchical analog where real space renormalization procedures are easier to define and to study: 
since its introduction for the classical ferromagnetic Ising model this framework has been used recently for 

many long-ranged disordered models, either classical like random fields Ising models [7^75] and spin-glasses |76l - l82j| . 
or quantum like Anderson localization models [ssl - l^ and quantum spin models [oll, . In this section, we consider 
the Dyson hierarchical version of the Long-Ranged Quantum Ising Chain 

H (23) 

i i<j 

where the properties of the long-ranged couplings Jqj, namely their decay with the distance (j — i) and their signs, 
have to be specified. 


A. Dyson hierarchical version of the Long-Ranged Quantum Ising Chain 

The Dyson quantum Ising Hamiltonian for A = 2" spins is defined as a sum over the generations A: = 0,1,n — 1 

n—1 




2 ") 


(24) 


fc =0 


The Hamiltonian of lowest generation k = 0 contains the transverse fields hi and the lowest order couplings ■ 
associated to the length Lq = 2° = 1 


H, 


(fc=0) _ 

( 1 . 2 " 




(0) X X 


(25) 


The Hamiltonian of next generation k = 1 contains couplings associated to the length Li = 2^ = 2 
2 " 


H)': = - ^ [E-3.4*-l'^4i-3^4*-l + "^1-3.4^41-3^4* + A^-2Ai-l<^-2<^-l + Jli-2^^^l^-2<i\ (26) 


r(fc=l) _ 
'(1.2") - 


r(i) 


r(i) 


and so on up to the last generation k = n — 1 associated to the length L„-i = 2" ^ couples all pairs of 

spins between the two halves of the system 


- E E If 


j * j 


i=l j=2"-i + l 

To mimic a power-law decay of the non-hierarchical model of Eq. 


(27) 


(28) 


one chooses the same dependence for the scale of the Dyson coupling of generation /c as a function of the length 

Lfc = 2'= 


j(fc) oc E = 2 

LI 


— ka 


(29) 


Within a given generation k, the couplings may be chosen uniform =2 ''“as in the Dyson quantum 

ferromagnetic Ising model with uniform or random transverse fields discussed in m or may be taken as random 
variables as in the Dyson quantum Spin-Glass discussed in 92|. In the following, we extend the real-space RG 
procedure introduced to construct the ground state of these models [U, 113 in order to construct also the excited 
states. 
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B. Diagonalization of the lowest generation k — 0 


The elementary renormalization step concerns the block two-spin Hamiltonian of generation k = 0 


— ^2i—lO'2i_i ^2iO'2j ^ 21 — 1 , 21^21 — 

The diagonalization within the symmetric sector 

H{2i-l,2i) I + + > = —{h2i-l + h2i)\ > —<72l*-l,2il-> 




H(2i-l,2i)\ -> — “>^2^-1,211 + + > +(^2i-l + ^2i)l-> 


yields the two eigenvalues 


=±\|{Atl, 2 ^? + i.h 2 ,-^ + h 2 .y 


corresponding to the eigenvectors 


l^ 2 i > =cos6»|'J-b-f >-bsin6»*J--> 

|Af+> =-sin6»|j-F-h >-bcosdfj - - > 


in terms of the angle 0 ^^ defined by 


cos(6»fi) = 


1-b 


h 2 i-l-\-h 2 i 


^/i4tl,2.r + (h 2 i-l+h 2 iy 


sin(dfj = sgn(J^°l^ 2i)\ 


1 - 


^2i- 1 +^2i 


4(4ti.2ir+ih2,-i+h2ir 


Similarly the diagonalization within the antisymmetric sector 

H(2i-l,2i) \ H-> = —(^2j-1 — h2i)\ -> —^2i-l,2j-> 

H(2i-l,2i)\ -1" > = ~J2i-l,2i\ -> +(^2i-l — h2i)\ -h > 

leads to the two eigenvalues 


( 30 ) 


(31) 


(32) 


(33) 


(34) 


(35) 


= ±\/(4ill.2i)^ + ih2^-l - h2iY 

with the corresponding eigenvectors 

\^t~ > = cosd^il -b - > Tsind^J - -b > 
1 X ^4 > = -sin0^j| -b - > -bcos0^j| - -b > 


in terms of the angle 9a defined by 


cos{0^) 


sin(0^i) 



^ _j_ h 2 i-i—h 2 i _ 

4i4f-l,2iy + ih 2 i-l-h 2 iy 
2 


= sgn(ji°li^2i)\ 


_ h 2 i-i—h 2 i _ 

4(4tl.2.)^ + (h2i-l-h2.y 

2 


(36) 


(37) 


(38) 


To construct the ground states, one chooses to keep only the two negative energy levels Af) and A^j [9l|,[93- Here 
we wish to keep also the two positive energy levels A^j^ and X^^ ■ As in the previous section, it is thus convenient 
to introduce a pseudo-spin T^( 2 i) ~ ^ choice between the subspace of positive or negative eigenvalues, and 
























another spin cr/j( 2 i) for the subsequent choice between the symmetric and the anti-symmetric sectors. More precisely, 
we relabel the four eigenstates as follows 

l^2i > = Vr{ 2 i) = = + > 

\^ 2 i > = l''■fl(2^) = +^<^R( 2 i) = ~ > 

\^ 2 i^ > = l''■fl(2^) = ~^^R{ 2 i) = + > 

l"^2i^ > — \'^R(2i) — ~^^R{2i) — ~ ^ (39) 

The block two-spin Hamiltonian of Eq. 1301 then reads 


H(2i-l,2t) — -^R(2i)TR(2i) - hR{2i)(^R(2i) 

where the renormalized field for 


(40) 


\J ('^2j-l,2i)^ + (^2i-l + h2iY + \J + (^2i-l - h2iY 

^R{ 2 i) = -^- 

is bigger than the renormalized field for 


(41) 


hR( 2 i) 


\J ('^2i-l,2i)^ + (^2i-l + h2tY — \J (</2i-l,2j)^ + (^2i-l — /i2i)^ 

2 

2/l2i-1^2i 

\l ('^2i-l,2i)^ + (^2i-l + h2iY + \J ('/2i-l,2j)^ + (^2j-l — /i2i)^ 


(42) 


In the following, we thus make the approximation that the pseudo-spins T^( 2 i) T “ 2"“^ emergent local 

conserved operators. As in the previous section, this approximation is possible only when the initial transverse fields 
hi and/or the initial couplings j'f'^ are distributed with continuous probability distributions, so that there isn’t any 
exact degeneracy between the fields flii( 2 i) of Eq. [TT] 


C. Renormalization of the couplings of higher generations k > 0 

We need to compute the initial operators erf in the new basis of Eq. [39] Within the approximation that the pseudo¬ 
spins T^( 2 i) emergent local conserved operators, one neglects all terms involving the flip operators {T^(^ 2 i) ’'''R( 2 i))’ 
and one obtains 


in terms of the constants 


^2i-l — ^2i-lT^(2i)>^R{2i) 

<^ 2 i — C2iCrfi(2j) 


(43) 


C 2 i-i = {cose^,smef^ + sine^icos 0 ^) = sgn( 






C 2 , = (cos 6»D cos e^i - sin 6»D sin 9^) = \^ 


1-f 






(44) 


As a consequence, the renormalized coupling between the spins <jR( 2 i) and cr/j( 2 j) is given by the following linear 
combination of the four initial couplings of generation k associated to the positions (2i — 1, 2 i) and (2j — 1, 2j) 


T t(^) I 2 I z I z z //1C^ 

j R(2i),R(2j) — ^2i,2j'^2iC2j + 'Hj(2i)'^2i-l,2:/'^2i-lC2j + TR(2j)’J2i,2j-l‘^2iC2j-l + TR{2i)^R{2j)’^2i-l,2j-l'^2i-lC2j-l[^0) 


r(^) 


(fc) 


t(^) 


So here in contrast with the short-ranged RG rule of Eq. [T31 the choice of the pseudo-spins values (T^( 2 i) = =*=1^ 
±1) will produce different amplitudes for the renormalized coupling in the different sectors. 
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D. Renormalization procedure for the non-hierarchical long-ranged model 

The block renormalization procedure explained above can be applied to the non-hierarchical long-ranged model 
of Eq. [531 The only difference is that in the RG rule of Eq. USj the initial couplings have not exactly the same 
properties: J 2 i, 2 j and J 2 i-i, 2 j-i are associated to the distance (2j — 2 i), whereas J 2 i, 2 j-i is associated to the sligthly 
smaller distance (2j — 2i — 1) and J 2 i-i, 2 j is associated to the sligthly bigger distance (2j — 2i -\- 1). Whenever the 
long-range nature dominates over the short-range case, one expect that the initial power-law model of Eq. |53| and 
its Dyson hierarchical analog have the same scaling properties at large scale. However the RG flows for the Dyson 
analog are usually easier to study analytically and clearer to analyze numerically. 


E. Iteration of the renormalization process 


When the renormalization process is iterated as explained in section III PI with the appropriate renormalization rules 
of Eqs im [4^ and |45l it should be stressed that not only the pseudo-spins of higher generation depend on the lower 
pseudo-spins, but also the amplitudes of the renormalized fields and couplings. More precisely, we may write the 
Hamiltonian after the first RG step 


H, 


( 1 . 2 " 


i=l 




R{2i)TR(2i) 


i=l 


T'R{2i)^R{2i) 


i<j 


R(2j) {T^R(2i) > ^R(2j) )^R{2i) ^R(2j) 


(46) 


after the second RG step 


^(1,2^) — ^ ^R{2i)'''R{2i) ^ ^^J^(4i)(''‘fl(4i-2))'^fi(4i))''’fl2(4i) (47) 

i=l i=l 

^n-2 

E l / 2 z \ z 't ( Z Z Z Z Z Z \ X cc 

IT-R? (4i) KT'R(Ai-2) j '^R{Ai) R{Ai) ~ -JR{Ai) ,iJ(4i) (.T'fl(4i-2)) '^R{Ai )) '^R{Aj-2 )) '’’fl(4j)) '^R'^ (4i)) '^R? ( 44 ) R{Ai) ^R(ij) 


i<j 


and so on up the last RG step 

2 " 


77(^,2") - “ X! “ X! ^fi''(4i)( 


'^R(4i-2) J '’'fl(4i))''’fi2(4i) 


i=l 


E v—V ( z z z z z z \ z 

“fl3(8q(T^2(gi_4), T^ 2 ( 8 j); 

1=1 

—H/j..(2n)({r^fc})T4j„(2") ~ (48) 

To make the link with Eq. |TJ one needs to expand the renormalized fields ^Irp and /i_r»( 2 ") in terms of the pseudo¬ 
spins of lower generations using the principle of Eq. [18] : each renormalized field of the second generation can be 
expanded into 2^ = 4 terms 


Q / 2 2 \ _ q{0i0) _i_ r)(lT0) z I q(0,1) z , o(l-l) 

^^R^{4iyR{4i-2) ^ ^R‘^ {4iy R{4i) ^^R‘^{4iyR{4i-2yR{4i) 


(49) 


each renormalized field f^fi3(8i)(Tfl2(8i_4), r^2(8q; r^( 4 j_g), T^( 8 i_ 4 ), t^( 8 j_ 2 )>'^^( 8 *)) generation can be ex¬ 

panded into 2® terms and so on. This rewriting will thus generate all possible products containing an arbitrary number 
of pseudo-spins as in Eq. [Tj 

As a consequence of the random sum structure of Eq. [48] over all blocks and all scales, one expects that at 
the middle of the spectrum near zero energy, two energy levels that happen to be consecutive (i.e. that have an 
exponentially small energy difference of order (N 2 2~^) with respect to the number = 2" spins) have completely 
different wave-functions labelled by completely different values of the pseudo-spins. 
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F. Properties of the eigenstates for the Dyson quantum spin-glass model 

To be more specific, let us now focus on the case of the Dyson quantum spin-glass model, where the initial couplings 
of generation k are random Gaussian variables of zero mean and of variance 

= = i (50) 

The parameter tr governing the decay with the distance has to be in the region ct > 1/2 in order to have an extensive 
energy. The initial transverse fields are taken uniform hi = h (but of course the renormalized transverse fields are 
random as a consequence of the RG rule of Eq. |35]). For an excited state, the random choice of the pseudo-spins 
T^( 2 j) = ± in the RG rule of Eq. |45] simply amounts to randomly change the signs of the couplings that are already 
of random sign. As a consequence, the RG flow is exactly the same as for the the ground state in another disorder 
realization. We can thus directly use the results concerning the RG flows for the ground state studied in as a 
function of the control parameter h : the constructed excited eigenstates are similarly either in the paramagnetric 
phase, in the spin-glass phase or critical with a finite dynamical exponent z{a). 


IV. CONCLUSION 

In this paper, we have proposed to construct the extensive number of local conserved operators characterizing a 
Fully Many-Body Localized quantum disordered system via some block real-space renormalization. The general idea 
is that each RG step diagonalizes the smallest remaining blocks and produces a conserved pseudo-spin for each block: 
for a chain of N spins, the first RG step produces ^ conserved operators associated to the blocks of L = 2 sites, the 
second RG step produces ^ conserved operators associated to the blocks of L = 4 sites, and so on. The final output 
for a chain of N spins is a hierarchical organization of the conserved operators with layers. We have explained 

why the system-size nature of the conserved operators of the top layers is necessary to describe the possible long- 
ranged spin-glass order of the excited eigenstates and the possible critical points between different FMBL phases. We 
have discussed the similarities and the differences with the Strong Disorder RSRG-X method [27l - l3lj| that constructs 
the whole set of the 2^ eigenstates via a binary tree of N layers. The block RG construction of the whole spectrum 
has been described for the random quantum Ising chain, first with short-ranged couplings and then with long-ranged 
couplings. In the Spin-Glass models, we have obtained that the excited eigenstates are exactly like ground states in 
another disorder realization, so that they can be either in the paramagnetic phase, in the spin-glass phase or critical. 
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